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Multivariate maxima of moving multivariate maxima 
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Abstract: We define a class of multivariate maxima of moving multivariate maxima, generalising the M4 pro- 
cesses. For these stationary multivariate time series we characterise the joint distribution of extremes and compute 
the multivariate extremal index. We derive the bivariate upper tail dependence coefficients and the extremal co- 
efficient of the new limiting multivariate extreme value distributions. 
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o 

p ( ; 1 Introduction 

' The M4 processes considered by Smith and Weissman (1996) constitute quite a large and flexible class of 

multivariate time series models exhibiting clustering of extremes. In the following generalisation we will 
introduce a copula in the array of innovations, rendering the new model more flexible than the previously 
p/| mentioned M4 processes, to which the new model reduces if the copula is the Frechet-Hoeffding upper 

p ' bound. The new family (M5) comprises multivariate time series models exhibiting clustering of extremes 

across variables and over time. We derive the limiting distribution of the vector of componentwise max- 
ima and the corresponding limit for the sequence of independent and identically distributed variables. 
Because of temporal dependence, these two limiting distributions are different, and the difference is quan- 
tified by the multivariate extremal index, which will be calculated too. 
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' Let {Z;.„ = {Zi_ 7lt i, . . . , Zi_ n d)}i>i.-oo<n<oo be an array of independent random vectors with stan- 

J> , dard Frechet margins and common copula C'z ■ A multivariate maxima of moving multivariate maxima 

■ (henceforth M5) process is defined by 

on ; 

Y nJ = max max a^ k ,j Zl,n-k,j , j = 1, • • • , d, n > 1, (1) 

, />1— oo<fc<-|-oo 

' where {a^k.j, I > 1, < k < oo, 1 < j < d} are nonnegative constants satisfying 



oo oo 



Y] Y] ai,k,j = 1, for j = l,...,d. 

xr 

H 

When Zi tn j = Z^ n , j = the M5 process is the M4 process considered by Smith and Weissman 

(1996). 

The common distribution Fy of Y n = (Y n> i, ■ • ■ , Y n _d) satisfies 



oo oo 

F Y ( 2/1 ,...,y d )=TT FJ F z (_?«_,..., , •(!.., 1 d, 

\a;,ft,l Cii, k ,d J 



Z — 1 fc — — oo 



and the relation for the corresponding copulas is 



oo oo 

c Y (Mi,...,^)=n n ^k w ,-,c m ), Uj €[o,i],j=i,...,d. 

l — l fe— — oo 
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This paper is concerned with the extreme value properties of the class of stationary processes defined 
in (p}. The M5 processes contribute to the modelling of variables moving together and generalise the 
important class of moving maxima discussed by several authors (Davis and Resnick, 1989; Deheuvels, 
1983; Hall et al, 2002; Smith and Weissman, 1996; Zhang and Smith, 2004; among others). This is the 
motivation of this paper which is organised as follows. 

By choosing a positive lower orthant dependent distribution function Cz in the domain of attraction 
of a max-stable copula C* (Joe, 1997) we will find a new class of limiting multivariate extreme value 
(MEV) distribution H for the vector M„ = (M„ ; i, . . . , M„ j£ z) of componentwise maxima from Yi, Y„. 

In the third section, we derive the multivariate extremal index of the M5 processes and illustrate the 
result with some choices of C* . 

Finally, we compare the bivariate upper tail dependence coefficients of Fy with the ones of the limiting 
MEV distribution H. 

2 Domains of max-attraction 

Let {Y„}„ > 1 be a sequence of independent random vectors associated to {Y„}„ > 1, that is such that 
Fy- = Fy, and M„ = ^M„.i, . . . ,M n ,cij be the corresponding vector of pointwise maxima. 

We shall assume that the copula Cz of the array of innovations in {Y„}„ > 1 is positive lower orthant 
dependent, that is, it satisfies the inequality 

d 

C z (ui,...,u d )>IJ« i) (u u -,u d ) e [0,l] d . (2) 

i=i 

Positive lower orthant dependence is implied by several dependence concepts (Joe, 1997, chapter 2)) and, 
in particular, is satisfied by associated random variables. The multivariate extreme value distributions 
are associated and therefore the following results can be applied when the array of innovations Z; „ has 
a MEV distribution. 

Following Smith and Weissman (1996), we present in this section the limiting distributions of the 
normalised vectors M„ and M„. 



Proposition 2.1 If C'z is positive lower orthant dependent and is in the domain of attraction of a max- 
stable C* , that is, 

C Z {U{ ,...,U d ' ) > C {Ut,...,Ud), 

n — ¥oo 

then Cy is in the domain of attraction of 

OO OO 

C( Ul ,...,u d ) = l[ [] C*K^\...,^ d ) (3) 

l—l k— — co 

Proof. We have 



OO OO 

/— 1 A;— -co 



and we can take the product of limits since 



oo 00 



lim J ' Cl(uT' k ' l/n ,...,u a d ^' n 



00 OO 



= ex P L 1 ™ E E lo &( C z(< 

l — l k— — oo \ l — l k— — oo 

and, by applying the Frechet lower bound in ([5]), for each k and I it holds 



1 • •) u d 



logCS(<' fc - l/,l ,...,^-" /n 



< 



n 



ai,k,j/n 



with X^ti I log I summable for each j. It then follows, from the Domminated Conver- 

gence Theorem, that 



00 00 



I™ n n ci(u?^ /n ,. 



ai,k,d/n 



i=l fe=-c 



00 00 



ex P E E log QimCg 



^ Z— 1 /e— — 00 



We have then, for Vr = (n, . . . , Td) € IR+, 



= C(iii, ...,Ud) . 



□ 



P M„,i <-,..., M, 

Tl 



><*< ^) =^y (e- Tl/ ",...,e- Td /") 



OO OO 



C*(e~ ri ,...,e~ Td ) = I ' C* (e- Tiai > h - 1 ,...,e- Tdai - k > d ) 

l — l k— — oo 



Let H denote the multivariate extreme value distribution with standard Frechet margins and copula 

OO OO 

C, i.e. H(xi, ...,Xd) = J C* ^e _:El lQi fc ' 1 , e _2: rf lQ, fc li ^ . We will now consider the corresponding 

l — l k — — OG 

limit iJ for the vector M„ = (M n i, . . . , M„ t d) of componentwise maxima from Yi, Y n . 



Proposition 2.2 // Cz *s a positive lower orthant dependent copula in the domain of attraction of a 
max-stable C* then 



lim P M nA < M„ td <—) = T[ C * 

7~1 7~rl ' 



1 = 1 



max ai^m - max a^ k ,dTd 

p — 00<&< + 00 p — 00</c< + 00 



Proof. It holds 

I ft ft 

lim P M„,i < M n . d < — 

n^oo V n T d 



(4) 



OO OO 



n c - « 1 



max ai^irx/n 



m<k<n—7n 



,...,e 



max a iiM T d /n N 

1 — m < A; < n — m 



l—l m— — 00 



OO OO 



™ cxp E E lo e Cz e 1 

m= — 00 \ 



max ai,k,xT\/n - max ai jk ,dTd/n^ 

-7n<k<.n—m „ l — m<k<n—m 

. .... c- 
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Now we can interchange the hrst two limits, by using analogous arguments to those used in Proposition 
2.1, since 



E 



m— — oo 



/ - max at^m/n 
log Cz I c ^ — 7n —^— n ~ g i 



max ai, k ,dTd/n 



oo 

* E 



log 



I d ^ 

~E, ^ a l,k,j r j/ n 
' •* l—m<k<n—m 

e i =1 



\ 



J 



oo d oo , 

El ^ — -s ^ — > d 

— > max oti k jTj < > — max max ai k -iTj 

n*— 1 l-m<k<n-m ' ' J J ^ n l<j<d l-m<k<n-m ' ' J J 



m— — oo j—1 



and this last sum is bounded for all n by X)fcL-oo ^ max { a i k jTj), which is summable in I (Smith and 

i<i<d 

Weissman, 1996). 
We have therefore 



lim P ( M„,i < M n . d < — 

n-too \ T\ Td 



exp( V lim V -logCde 1 



max OLifc\T\jn 



— m<k< 



max ai.k.dTd/n" 

l — m<k<n—m 



\l — 1 m— — oo 

and we have to prove that 

oo 



lim — log C 



n— foo * — ' 71 

m— — oo 



max a L k,\Ti/n - max oti,k,dT d /n 

'n | g 1— in<k<n — in g 1 — m<k<n — m 



max a;,fc,iri 
logC* e -°o<fe<+°° 



..,e 



max a;,fc,dT d 

— oo<fc<+oo 



For fixed /, let bi & j = ai k jTj, j = 1, d, and bi k = max bi & j (which are summable in k and I). 

j=l,...,d 

Suppose that bi_k.j is maximized when ft = k*(l,j), bi t k is maximized when k = k*{l) (not necessarily 
unique) and assume, without loss of generality, that k*(l, 1) < ... < k*(l,d). Break the above sum into 
three sums Sn\ 1 — 1,2,3, corresponding to 1 — min{/c*(Z, 1), &*(/)} to < n — max{/c*(Z, <i), &*(/)}, 
to < 1 - min{fc*(Z, 1), fc*(/)} and to > n — max{fc*(/, d), k*(l)}. 
For large n we have 



5 j. 



1 



n-max{fc* (l,d),k* (i)} 



E 



logCg e 



I „ l-m<fc< 



max ai t k,m/n 



l-m<k< 



max ai t h,dTd/n^ 



m=l-mm{fc*(i,l),fc*(Z)} \ 

n - max{F(i, d), fc*(Z)} + min{fc*(Z, 1), k*(l)} 



J g(^n ^g-ai,fc*(l,l)lTl/n e -ai,k*{l,d)dTd/n^j 



which converges to 



log C* (e _a, > fc *( i > 1 > iri e ~ Q i.fc* d,d)d T d^ 



= \ogC* 



max az,fc,iTi 

, — oo<fc<+oo 



...,e 



max a Lk ,dTd s 

-oo<A;<+oc 
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Otherwise, by applying @, we can write 



si 



< 



d 

1 E E 

n z — ' 4-* l- 



max m.kjTj 

71 L — ^ ^ — ' l — m<k<n—m 
m<l-fc*(/) J=l 



1 

n 



E E, 



l-Jb»(0<m<X-min{fc*(J,l),fc*(0} i =1 



max &l,k,jTj 

7n<k<n—m 



< 



max max oti,k,jTj 

71 * — ' 1— m</c<n — m ^ — 
m<l— fe* (i) 



d 
n 



E 



l-/:*(Z)<m<l-niin{fc*(Z,l),/c*(Z)} 



max max oti^jTj 

l-m<Kn-mj-l,...,tJ ' ' 



and 



< - max max a ( fc ,-r,- + - (min{/c*(J, 1), k*(l)} — k*(l)) b, wn 

71 A — ' l—m<k<n—mj — l,...,d 71 
m<l— k* (I) 



m>n— fc* (Z) j=l 



max oiikjTj 

m<k<n—m 



1 



E 



-max{fc* (l,d),k* (l)}<m<n — k* (I) j=l 



y~] max ai.k,]T] 



< 



d 
n 



Emax max a; t ,-t,' 



i>n—k* (I) 



l-m<Kn-mj-l,...,(f 



d 
n 



E 



i-max{fc* (Z,d),fc*(Z)}<m<n — 



max max Q£i,k,jTj 

1 — m<fc<n— m j—l,...,d 



< 



\^ max max Qti,k,j T j H — (max{A;*(Z, d), &*(/)} — fc*(Z)) bik*(l) 

^ — ' 1 — m<k<n— m j — l,...,d TL 



m>n—k* (I) 



The second terms in the above upper bounds for 



i = 1,2, tend to zero and the same holds for 



the first terms, by Lemma 3.2 in Smith and Weissman (1996). 



□ 

For the particular case of C*(u%, Ud) = min Uj, (ui, . . . , Ud) S [0, l] d , the above result agrees with 

l<j<d 

(3.5) in Smith and Weissman (1996), that is 



lim P I M n x < M n d < — ) = exp < - max max aik,jTj 

Tl T~dJ ^—^ -oo<fe<oo l<j<d J J 



If we choose C*(ui, Ud) — T[j=i u ji • • • j u d) € [0, l] d , we find again the result in case two of 
the last example in Martins and Ferreira (2005), 



M n l <—,..., M n d < — ) = exp < max ai,k,j T j 

In this last work it is considered only Cz = C* and the above two choices for C* are treated directly. 

By calculating the limits in ([3]) and for a given initial copula C* , the above propositions enable 
us to obtain a larger class of MEV models. 

Some choices for C* are available in the literature. For instance, Caperaa it et al. (2000) gives 
sufficient conditions for an Archimedean copula Cz to be in the domain of attraction of the Gumbel- 

Hougaard or logistic copula C*(u±, Ud) = exp ^— (Y^—i ( — l°g u j) Q ) ■ More examples of Cz in 

the domain of attraction of a MEV copula C* can be found in Demarta and McNeil (2005) and Hiisler 
and Reiss (1989). 

3 The multivariate extremal index of M5 process 

In this section we will extend the results about the multivariate extremal index of the M4 process, as a 
corollary of the propositions 2.1 and 2.2. 

We first recall the definition of the multivariate extremal index function 9(t) = 9{t\, . . . , Td), t G 
that relates the MEV distribution functions H and H and which was introduced by Nandagopalan (1990). 

A d— dimensional stationary sequence {Y n }„>i is said to have a multivariate extremal index 9(t) 6 
[0,1], r £ IR+, if for each r = (ri,...,Td) in Hi, there exists = 
{Un\\- ■ • , u { J d J), n>l, satisfying 



nP(Y lj >u™) — > - . j 1 ,/. 

J n— >oo 

P(M„ < uM) — > t(t) G (0,1] and P(M„ < u^) — > 7 (r) 6 



Proposition 3.1 // Cz is a positive lower orthant dependent copula in the domain of attraction of a 
max-stable C* then the multivariate extremal index of the M5 process {Y„} defined in (J^l is given by 



]Tio g c* 



max ai t k,m - max a^ k ,dT d 

-00<&< + OG p — OG<fc< + 00 



(t%, . . . , Td) — oo oo 



J2 J2 l °8 C * (e- ai - k - lT \... 1 e~ ai - k - dTd ) 



1=1 k= 



and the extremal index of {V n ,j}n>i is 

oo 

6j = V" max ai.k.j , j = l,...,d. 



-oo<fc<oo 

1=1 



(t) r- 

Proof. By applying Proposition 2.1 and Proposition 2.2 with u n 3 - = j = 1, d, we find 



oo oo 



7(r) = J ' C* (e- ai ' fc - lT \...,e- ai ' fc < d7 



Z=l fc=-c 

and 



SZ. f - max aj t,m - max a^ k . d T d \ 
7 (t) = 7(t) 8 W = [](7* e 7 ... ;6 -°o<fc<+oc, \ 

i=i ^ ' 

that leads to the results for 0(ti, . . . , t^) = ! os ?f T | and 0,- = lim #(ti, . . . , rA 

log 7^; r ._ > .o+,ie{i <i}-{j} 



□ 



Therefore the copulas C*, C and the copula C of the limiting MEV distribution H are related by a 
multivariate extremal index 0{t\, . . . , r^) through 

v 01 '°g"l '°8 "d \ OO OO „, log hi 

/o/ \ (r<( Vfi i/^nV Sl " d TT TT (n*( Q i,<=,i/ e i «i,m/ 9 aV 81 

C(ui,...,Ud) = ^(u/ ,...,u d ' )J =_[_[ _[_[ [C (uj ,-,?V )J 

/— 1 fc — — CJO 

The proposition 3.1 leads to the results of Smith and Weissman (1996) and Martins and Ferreira 
(2005) when C* is the copula of the minimum and product, respectively. 
Otherwise, if we take for instance the Logistic copula C* we find 

l/a 



1=1 \j=l 

0(Ti,...,T d ) 



EV* ( max ai,k,jT 3 



, i/« 

oo oo / a 

EE E ( a i,k,j T j) c 

1=1 k—oo \j— 1 



4 The tail dependence of M5 process 

For a random vector X = (Xi, . . . , Xd) with continuous margins F d and copula C, let the bivariate 

(upper) tail dependence coefficients parameters be defined by 



A 



\f} = timPWXi) > u\F f (X f ) >u),l<j< f < d. 

J J Ujl 



The tail dependence coefficient characterizes the dependence in the tail of a random pair {Xj , X.y ) , 
i.e., Xj ), > corresponds to tail dependence and X?), = means tail independence, and can be defined 
via the copula of the random vector which refers to their dependence structure independently of their 
marginal distributions. 

It holds 

Affi=2-lim tog ^' (u ' u) , 

JJ ufl logW 

where Cjji is the copula of the sub- vector (Xj,Xj/). 
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(C) (C) 

In this section we will relate A^; with A^-,, where C is copula of H and C is the copula of H. We 
first remark that, for each (uj,Uji) £ [0, l] 2 , we have 



where 0(rj,Tj>) is the bivariate extremal index of {(Y n j ,Y n ji)} n >i. This relation enables us to compare 
the tail dependence parameters X^ C J with Xj C ), through the function 0(rj,Tj'). 

Proposition 4.1 If C and C satisfy fSJ] then 
(a) xf)=2 + e{j-.-^)\ogC j , j .{e-^,e- 1 ^'), 



Proof. From the spectral measure representation (Resnick, 1987) of the copula C, with measure W, 
we get 



\f), = 2 - lim B " v ' ' = 2 - / max^.io^dW^wi,...,^) = 2 + logC^v (e^e" 1 ). (6) 



,(C)_ 9 logC , j, J -/(u,u) 

«ti bgu y 5d 

By using ((5|) and the homogeneity of order of the multivariate extremal index, it follows that 

A (C) = 2 _ lim \QgC, r (u,u) = 2 _ lime(Z log« ^logtx logC^^,^') 

JJ ufl logU ufl 0j By logW 

1 1 Js rf max { — tr^> — w - r ^)dW{w 1 ,...,w d ) 
= 2-0( — , — )lim— 3 - 

= 2-ed-,±-)f m^{^L,^}dW( Wll ...,w d ) 



= 2 + 0(1 J-)logC^(e- 1 ^ ! e- 1 /^) ) 

which has as a particular case. To obtain the second statement we combine the first with ^ and use 
the max-stability of Cj^i , as follows: 

xf) = 2 + log C jtj , {e-\ e- 1 ) - log C,. r {e~\ e^) + ±-) log C^e" 1 ^, e -V^) 
= A (C) , l0£ (^'( e ' 1/ ^ e " 1/9j '))^^ ) 

Qj' ( e_1 ! e_1 ) 

= Aff+log ^' (e ^— 1 

□ 



We will now apply these relations to the particular case of M5 processes, enhancing the effect of the 
multivariate extremal index in the tail dependence. 
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Proposition 4.2 Let {Y„}„ > 1 be a M5 process defined as in (QJ) and such that Cz a, positive lower 
orthant dependent copula in the domain of attraction of C* . Then, for any 1 < j < j' < d it holds 



oo oo 



l — l k— — oo 



oo oo 



(b) aS2= 2 + £ lo g q,,(e-^^/^,e-«^.W^ 

2=1 k— — oo 

Proof. By using © and the Proposition 2.1, we obtain (a); by using the Proposition (4.1)-(a) and 
Proposition 2.1, we get (b). □ 



If C*(ui, Ud) = min Uj, (u-i, Ud) S [0, l] d , we obtain 

l<j<d 



oo oo 



A L-' 2_ E E max { a ',fc,j: a i,fcj'} ( 7 ) 



2=1 /c— — oo 



in general greater than zero, which agrees with the result (2.10) presented in Heffernan et al. (2007). 
We can also say from (b) that, for this copula, 



OO s 

Xy/, = 2 — V* max \ max max a Lk 

■>'■> ^ -oc<fc< + oo J J -oo<fe<+oo 

2 = 1 



and it holds Xj2 > Xjjl if and only if 



OO OO OO s *. 

maxjaj.fc.j, ai. k .f} - Y^max <^ max ai k j/6j, max a>i t j> /6j> \ > 0, 

* — * * — * ^ — * — oo</c<+oo — oo<fc<+oo 

2=1 fc=-oo 2=1 k ~ ~ ~ ~ ^ 



where 6j = N max cti,k,j, j = 1 5 • • 



-oo<fc<oo 
2=1 

(C) (C) (C) (C) 

We can easily construct examples of M4 processes for which Xj ■; > XjJ or Xj ■; < AV;. 

The results for this class of moving multivariate maxima show us that if we obtain or estimate a tail 
dependence parameter of the common distribution of the variables in a stationary sequence we don't have 
necessarily the corresponding parameter in the limiting MEV model C . The multivariate extremal index 
of the stationary sequence can increase or decrease the tail dependence and the extremal coefficients of 
the limiting MEV model C arising from the i.i.d. sequence. 

Since \\ -', = AV, , the result in ([7]) says that in the M4 processes the variables Y n j, j = l,...,d, 
are in general asymptotically dependent. For the M5 processes we can choose C* in order to produce 
variables asymptotically independent. Take for instance C* such that Aj ; •, = and ai.k.j = czi,k,j'- We 

then find xQ = xfP = 0. 

The Proposition 4.2. also points out that even for a choice of a copula C* with symmetric tail 
dependencies, the values of the signatures a^k.j can lead to asymmetric tail dependencies in both copulas 
C and C, that is, to different values of Xjj/ for different choices of (j, j'). 



The results in the above proposition are translations of the classical result A = 2 — e for Cjji and Cj ji 
(Nelsen, 2006) In fact, if we define the extremal coefficient of the MEV copula C as the constant ec such 
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that C(u, ...,u) = u ec for all u S [0, 1], then from 

nt \ ( at lo % Ul l °S u d, f -logUjW ; » \ 
C(ui,...,Ud) = exp -8{— — — - — ) / max dW{w 1 , ...,w d ) , 

\ Ol Vd Js d 1 <3< d Vj J 

we find 

ec = -e(±...,±)logC(e-^,..., e -V<>«) 

and, in particular, 

e (? = -log(7(e- 1 ,...,e- 1 ). 

(C) 

The relation between Cjj* and C*y enables now to obtain the Proposition 4.2. from A^- -; = 2 — 
and Ag) = 2 - ec w , 

Other features of the M5 processes such as a directory of tail dependence coefficients for different C* 
and signatures (Xi,k,ji illustrating the range of dependence structures, the model selection and estimation 
are key directions in a future research. 
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